On electron (anti) localization in graphene 
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We discuss localization properties of the Dirac-like electronic states in monolayers of graphite. In 
the framework of a general disorder model, we identify the conditions under which such standard 
localization effects as, e.g., the logarithmic temperature-dependent conductivity correction appear 
to be strongly suppressed, as compared to the case of a two-dimensional electron gas with parabolic 
dispersion, in agreement with recent experimental observations. 



After several decades of predominantly application- 
driven studies, graphene has finally been recognized 
as a unique example of the system of two-dimensional 
fermions with a linear dispersion and pseudo-relativistic 
kinematical properties. The recent advances in microfab- 
rication of graphitic samples that are only a few carbon 
layers thicki^ have made it possible to test the early the- 
oretical predictions of the anomalous properties of this 
system^. 

The most striking experimental observation up to 
date was that of an anomalous quantization of the Hall 
conductivity^^ which is characteristic of the pseudo- 
relativistic nature of the quasiparticle excitations in this 
system. The other properties manifesting this peculiar 
single-particle kinematics have been revealed by mag- 
netotransport measurements, including the \/~B depen- 
dence of the energies of the (non-equidistant) Landau lev- 
els and the intrinsic 7r-shift of the phase of the Shubnikov- 
de Haas oscillations^*^. 

In a (nearly) degenerate semimetal such as graphite, 
the Coulomb interactions are expected to play an impor- 
tant role, largely due to their poor screening^. Besides, 
any interplay between the Coulomb interactions and dis- 
order is likely to further modify the behavior of the ide- 
alized (clean and non-interacting) Dirac fermion system. 
Thus far, however, experiment has not yet provided a 
direct evidence of any interaction-induced phenomena, 
some of which were predicted in the recent years&i. 

Also, contrary to the situation in the conventional 
2DEG, the data of Refsi^ do not seem to manifest any 
pronounced weak-localization effects, either. Motivated 
by these observations, in the present work we set out to 
study electron localization in graphene. Although this 
topic has already attracted some recent attention^, we 
shall demonstrate that the analysis of Ref£ is neither 
complete, nor (curiously enough, considering that the as- 
sertions made in Ref£ are mostly correct) do the argu- 
ments presented in that work appear to be technically 
sound. 

The electronic band structure of graphene is charac- 
terized by the presence of a pair of inequivalent nodal 
points at the wave vectors Ki,2 = 27r/9a(±Vo, 3) where 
a is the lattice spacing. At these two and the four other 
points in the Brillouin zone obtainable from K 1,2 with 
a shift by one of the reciprocal lattice vectors Qi,2 — 
27r/3a(v / 3, ±1), the valence and conduction bands touch 
upon each other as a pair of opposing cones with the 



opening angle given by the Fermi velocity vf- 

In the leading approximation, quasiparticle excitations 
in the vicinity of the nodal points (hereafter referred to 
as valleys) can be described by the Dirac Hamiltonian 3 



H = v F [l 



ffyky] 



(1) 



where cr, is the triplet of the Pauli matrices acting in the 
space of spinors = (ua,Vb) composed of the values of 
the electron wave function on the A and B sublattices of 
the hexagonal lattice of graphene, whereas the triplet Tj 
acts in the valley subspace. In the absence of magnetic 
field, the Hamiltonian (1) remains a unity matrix in the 
physical spin subspace. 

In the presence of disorder, the quasiparticles experi- 
ence both intra- and inter-valley scattering. Upon aver- 
aging over disorder, the most general form of the elastic 
four-fermion vertex induced by disorder can be described 
by the expression 
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where we used the customary "g-oloqical" notations 
to denote the processes of forward scattering between 
fermions from the same (54) and different (52) valleys, 
as well as those of "backward" scattering (<?i) with the 
momentum transfer close to K\ — K2. Besides, we have 
included the possibility of "umklapp" scattering (g 3 ), in 
which case the total momentum of two fermions changes 
by 2(^1 — ^2) = K2 — K1+Q. A justification for this (not 
immediately obvious, considering that the total momen- 
tum changes by only a fraction of the reciprocal latice 
vector) extension of the disorder model will be discussed 
below. 

In the framework of the self-consistent Born approxi- 
mation (SCBA), the effect of disorder on the quasipar- 
ticle spectrum is described by a self-energy which obeys 
the equation 



dq 

(2^y 



■TrWG R (u,k + q) (3) 



where the retarded Green function is given by the expres- 
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ajj, a _ - /i + 1 + g g^fcx + t z g o-^fc^) 
G K/Cj " (c- M + E«) 2 - V 2 fc 2 

that includes a chemical potential /i which allows one to 
account for a variable electron density. 

Notably, the solution to Eq. (3) turns out to be diagonal 
in the valley subspace, £ R (0,0) = 171 ® 1, and propor- 
tional to the inverse elastic lifetime 
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HVp 



(.94 + .9i) 



(5) 



where vp is the density of states at the Fermi energy. 
In the low-doping limit (// < 7), a finite value of vp = 
— J dk I (47r 3 )7mG(/i, k) is dominated by disorder. A sim- 
ilar phenomenon has been extensively studied in the con- 
text of normal quasiparticle transport in dirty c£-wave su- 
perconductors where /i = and the spectrum possesses 
an exact particle-hole symmetry^. In both cases one ob- 
tains the self-consistent DOS as vp — (7/271") \n(vp/aj). 

However, most of the data of Refs. 1 ^ pertain to the 
regime of relatively high dopings (// >> 7), as indicated 
by, e.g., the measured mean free path, which was found to 
be of order ~ Ifxm at electron densities n e ~ 10 13 cm~ 2 . 
In this regime, the value of vp = kp / (2-kvf) is controlled 
by the finite radius of the Fermi surface kp = (i/vp — 
(47rn e ) 1 / 2 and is only weakly affected by disorder. 

Therefore, the double-pole Green function (4) given by 
a four-by-four matrix can be well approximated by a pair 
of two-by-two matrices containing single-poles 



Gi, 2 (w,A;) 



a x cos 



b k ± by sin <f> k 



2(w - & + 17) 
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where £/. = vpk — [i and <pk = tan -1 k y /k x . 

Even more importantly, the presence of a large param- 
eter /z/7 facilitates a systematic account of quantum in- 
terference corrections to the SCBA results. We mention, 
in passing, that in the low-doping (/1 < 7) limit the only 
parameter that can (at least, in principle) be used for the 
analysis of such corrections is the (inverse) number of val- 
leys, whose actual value is, of course, N v = 2. The same 
caveat plagues the analysis of dirty d-wave superconduc- 
tors where the strength of the conductivity corrections 
is governed by the number of pairs of opposite nodes of 
the order parameter, whose physical value is again equal 
two. 

Luckily, in the high-density regime that, so far, has 
been probed in graphcne, the leading quantum correction 
to the zeroth-order (Drude) conductivity stem from the 
standard single-Cooperon (fan-shaped) diagram. From 
the technical standpoint, however, the calculation of the 
corresponding correction appears to be somewhat more 
involved due to the matrix structure of the Green func- 
tion (6) and the vertex (2). 

The Drude conductivity itself is given by the stan- 
dard expression o§ — (e 2 /h)T[i/j where the renormal- 
ization factor T accounts for a ladder series of vertex 



corrections associated with one of the two current op- 
erators J = up (cos <pk, sixKpk) inserted into the fermion 
loop in the diagrammatic representation of the Kubo for- 
mula. Being given by a non-singular (angular momentum 
m = 1) diffusion mode in the expansion over the angular 
harmonics e™^ fc , T is a function of the parameters gi. 

Expanding the expression for the Cooperon in the same 
basis as that used in Eq.(2), we obtain equations for the 
corresponding amplitudes Ci, each of which is a matrix 
in the direct product of two bi subspaces 

Cx = g x + g 2 H 12 Ci + giH 21 C 2 
C 2 — g 2 + g 2 H 12 C 2 + g\H 2 \C\ 
C3 = .93 + g^H\\Cz + g-iH 22 C4 
C 4 = g 4 + giHud + g?,H 22 Cz (7) 

where JTy(w,$) = / j0pG?(e + u>/2,k + q/2)Gf(w/2- 
e, q/2 — k) denotes a convolution of a pair of Green func- 
tions. Computing the latter at uj,VFq << 7, we obtain 

#11,12 (w,g) = #22,21 (w, 9 = — ^ 1 + — - 1®1 

47 27 167^ 

1 iuj v 2 F q 2 3cos 2 </> fc +sin 2 ^ fc 

-2 (1 + 2^ " W 4 °* 

T±d + £ - ^os 2 fc +3sin 2 ^ ® *,](8) 

It can be readily seen that Eqs.(7) split onto two pairs 
which only couple C\. 2 and C34, respectively. Their so- 
lutions read 



Ci >2 (w, 9) 



27d j2 1 ® 1 — b x (g) b x + b y ® b y + b z ® b z 
kvf {94 - g 2 )/(g 2 + 51) + v 2 F q 2 /16j 2 - iuj/2j 



- 2703,4 1 ® 1 - <Jx ® cr x - a y (g) <j y - a z ® a z 

C 3 ,4(w,9) = 7 T77 x , 2 2 , r a 2 ( 9 ) 

™f (31 - .93)/(53 + 34) + 149 /167 - Jw/27 

where a = (gj + gig±)/{g4, - .92 + 2gi)(gi + g 2 ), c 2 = 
(5254 + 5152 ~ g 2 + gf)/(g4 -92 + 2g\){g\ + g 2 ), c 3 = 
(.935i+53.94)/(5 , i+.93)(53+.9 4 ), and c 4 = (gm+g 2 )/ (.9i + 
53)(33 + 54)- 

The quantum conductivity correction (including spin) 
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■Tr[C 1 (0,q) + C 4 (0,q)] (10) 



involves the C\ and C4 components of the Cooperon, 
whose contributions turn out to be negative and positive, 
respectively. The logarithmic temperature-dependent 
part of Eq.(10) can be cast in the form 



(Jo - = — — In 



2e 2 max[r /7, \g 4 - g 2 \/(gi + g 2 )} Cl 



irh maxlT^/j, \gi - g 3 \/(gA + gz)] Ci 



(11) 



where we introduced an inelastic phase relaxation rate 
T$(T) which provides a cutoff in the momentum integra- 
tion and diverges at T — > 0. 
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The analysis of Eq.(ll) reveals that, in the absence of 
a fine tuning between the amplitudes of backward scat- 
tering and umklapp processes (gi ^ (73) the C4 Cooperon 
always acquires a gap ~ 7I91 — 93 1/(94 +93)- On the other 
hand, the C\ mode remains gapless, provided that all the 
forward scattering processes are controlled by the same 
amplitude (i.e., g 2 = 94). 

Conversely, making the C\ mode gapful and inverting 
the sign of the conductivity correction would only be pos- 
sible under the condition \g\ — g%\ << \g<i — g$\ which is 
unlikely to be satisfied for any realistic impurity poten- 
tial that yields equal amplitudes of the processes of intra- 
and inter-valley forward scattering. 

The antilocalizing behavior predicted in Rcf. 8 for 
9i « 92 = 94 and g 3 = (ci = 1/2, c 4 = 1) 
can only occur at intermediate temperatures (namely, at 
79i/92 < r^(T) < 7), whereas at still lower tempera- 
tures the overall sign of (11) reverts to negative, thereby 
suggesting the onset of rather conventional weak local- 
ization. 

The above conclusions apply to the general disorder 
model (2). However, in the situation where disorder is 
realized as a random distribution of impurities with a 
concentration n, and a (short-range) potential u(q), it 
suffices to introduce only two independent parameters 



92 = 94 



n lU 2 (0) 



1 + (tw f m(0)) 2 



9i = 93 



mu 2 ^ - K 2 ) 
1 + (txv f u(K x - K 2 )f 



(12) 



The above expressions represent a T-matrix computed 
for an arbitrary strength of disorder, the customary Born 
and unitarity (where the scattering phase approaches 
7r/2) limits corresponding to u — > and 00, respectively. 
[In the case of a genuine long-range (unscreened) impu- 
rity potential, the dependence u(q) on the transferred 
momentum makes the explicit formulas for gi more in- 
volved, though.] 

Provided that the relations (12) between the param- 
eters gi hold, one obtains c\ = C4 = 1/2, and the log- 
arithmic term in (11) vanishes as a result of the exact 
cancellation between the contributions of the localizing 
(Ci) and antilocalizing (C4) Cooperon modes. It has to 
be stressed, however, that such a strong suppression of 
(anti)localization would only be possible due to an open- 
ing of the umklapp channel. In turn, the latter requires 
an emergence of a crystal superstructure with the wave 
vector 2{K 1 - K 2 ) = (2/3)(Qi + Q 2 ) (or equivalent). 

While an isolated sheet of weakly-interacting graphene 
would apparently lack such a superstructure, it is con- 
ceivable that the latter might emerge if a commensurate 
substrate were used during the process of microfabfica- 
tion. Besides, a commensurate corrugation could occu»i£ 
due to the Coulomb correlations that can induce spatially 
periodic patterns of the electron density itself. The pos- 
sibility of a spontaneous formation of such charge density 



wave states has long been discussed in the general context 
of degenerate semimetals and, specifically, in grapheneii. 

Next, we comment on the technical details presented 
in Refi& where the first prediction of the antilocalization 
behavior for 172 = 94 and g\ = g% — was made. In 
essence, instead of working directly with the matrices in 
the sublattice subspace, the authors of Ref£ used single- 
component (projected) Green functions G(uj, k) = 1/ {u> — 
£fc + i"f) complemented by matrix elements of a (purely 
forward-scattering) impurity potential between different 
Bloch wave functions. 

To that end, in Ref. 8 the latter were chosen in a par- 
ticular (asymmetrical) gauge 



^ lfc =-L(±l, 0, 0), 

fhk = ^=(0, 0, e**, ±1) 



(13) 



The aforementioned matrix element of the impurity po- 
tential is then given by a complex-valued expression 
<p\u\k >oc (1 -fe'^-'W). 

Next, the authors of ReffS asserted that the equations 
for the Cooperon mode involve the disorder-induced ver- 
tex 



W k -k,p- P — rii < p\u\k >< -p\u\-k > 
oc e^-^l + cos(0 fc - <f> p )} 



(14) 



which describes scattering from the Bloch states with 
momenta k + q/2 and —k + q/2 into p + q/2 and —p + q/2 
in the limit q — * 0. 

It was then argued in Rcf. 8 that the full Cooperon am- 
plitude inherits the phase factor from Eq.(14) and, there- 
fore, becomes negative for k ~ — p where e l ^ fc ~^ p ) ~ — 1, 
which configuration of the momenta provides a dominant 
contribution to the weak localization correction. Thus 
obtained, the negative sign of the Cooperon amplitude 
was claimed to be instrumental for the onset of antilo- 
calizing behavior. 

The validity of the above argument can be disputed 
by observing that the choice of a different (symmetrical) 
gauge for the Bloch wave functions 

V4 fc = -^(±e-*/ 2 , e^/ 2 , 0, 0), 

^ = 4|(0, 0, e^/ 2 , ±e"^/ 2 ) (15) 

would result in a purely real matrix element < p|u|fc >oc 
2 cos ^ fc ~^ p ; the use of which yields Eq.(14) without the 
said phase factor, hence the (deemed to be crucially im- 
portant) sign change of the Cooperon amplitude does not 
seem to occur. 

Although the observed gauge dependence of the 
Cooperon is perfectly consistent with its being a gauge 
non-invariant two-particle amplitude, the conductivity 
correction is, of course, supposed to be gauge invariant. 
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A resolution of such apparent contradiction goes as fol- 
lows. 

In fact, the vertex that ought to be used in the 
construction of the Cooperon is the particle-hole (ex- 
change) amplitude W k , P , P ,k which is manifestly gauge in- 
variant and given by Eq.(14) without the phase factor 
in question. In this regard, it is important to realize 
that the amplitudes Wk,-k,p,-p = W qtq - k . P , q - P \ q ->o and 
W k . p,p,fc|fc— >-p represent two non-commuting (and, in this 
particular case, unequal) limits of the general vertex (2). 

On a side note, it might also be tempting to try to 
explain the experimentally found suppression of localiza- 
tion by the presence of the factor [1 + cos(<fik — <P P )] van- 
ishing for k = —p in Eq.(14). Notice, however, that the 
Cooperon built out of such a vertex lacks this factor, since 
the gapless pole develops for only one (namely, m = 0) 
angular harmonic, whereas the other two harmonics that 
constitute the above angular-dependent factor (m = ±1) 
receive only a finite enhancement by a factor of two. 

Thus, an obvious suppression of the backward scatter- 
ing off of an individual impurity (which fact has often 
been mentioned in the context of electron transport in 
ID carbon nanotubes) does not necessarily imply the ab- 
sence of (anti)localization in the deep diffusive regime in 
2D graphene. 

The true origin of the antilocalizing behavior exhib- 
ited by Eq.(ll) for |<?i — g 3 \ « \g 2 — g±\, can be traced 
back to the negative signature of the expansion of the 
Cooperon mode C4 over the product basis a a <S><Jb, as op- 
posed to that of C\ (see Eqs.(9)). Interpreting the sublat- 
tice index as a fictitious spin one-half (Si ,2 = 0X2/2), one 
can associate the C4 component with the singlet mode 
C*4 oc (1 ® 1 — |(Si + S2) 2 ) which, in the spin-orbit con- 
text, is known to be a common source of antilocalization. 

In this regard, the present anilocalizing behavior is 
similar to that predicted in the case of 2DEG with a 
strong spin-orbit coupling of either Rashba or Dressel- 
haus kincU£ 

H so = a(a x k y - a y k x ) + (3(& x k x - a y k y ) (16) 

The details of the calculations presented above differ from 
those pertaining to the case of the spin-orbit coupling (for 
one thing, except for a splitting, the linear in momen- 
tum SO-coupling does not alter the fermion dispersion, 
whereas in the case of graphene it is the linear term which 
is solely responsible for the dispersion). 

Yet another argument invoked by the authors of Rcf. 8 
to support their prediction of the antilocalizing behav- 
ior intermediate temperatures exploits the notion of the 
Berry phase associated with the electron wavefunctions 
(13). 

Indeed, the first Chern number (" vorticity" ) 

*=<y <^jkiv fc iv* >, (17) 

where the momentum space integral is taken along any 
contour enclosing one of the Dirac points if 1,2, yields 



$1.2 = — 7r for both wavefunctions (13). In Rcf. 8 (sec 
alsoi^) this observation was used to assert that the inter- 
ference between any pair of counter-propagating electron 
waves becomes destructive, thereby giving way to the an- 
tilocalizing behavior. 

However, the Berry phase (17) would be identically 
zero, if the wave functions (14) were used instead of (13), 
which suggests that its relationship with the Cooperon is 
somewhat more subtle. To that end, it is worth mention- 
ing that the two Dirac species emerging in the contin- 
uous (low-energy) description of the original lattice sys- 
tem must have opposite chiralities^*, as required for com- 
pliance with the ubiquitous Nielsen-Ninomiya ("fermion 
doubling" ) theorem^. In turn, this observation seems to 
suggest that, even if the Berry phases of the individual 
species happen to be finite ($1,2 = ±tt), they should still 
add up to zero. 

A cancellation of the overall Berry phase would also be 
consistent with the vanishing of the quantum conductiv- 
ity correction in the presence of equally strong intervalley 
backward and umklapp scattering, which, in effect, make 
the two Dirac points indistinguishable (the shortest re- 
ciprocal lattice vector becomes equivalent to K\ — K<i). 

Recently, the role of the Berry phase has also been 
brought up in the context of the spin Hall effect gen- 
erated by a spin-orbit coupling. Interestingly enough, 
the dependence of the spin Hall conductivity upon the 
coefficients a and j3 from Eq.(16) (and, especially, its ex- 
act vanishing at a = ±/3) can be established with the 
use of the Bloch wave functions in either asymmetrical 
(similar to Eq.(13)) or symmetrical (similar to Eq.(14)) 
gauge, where the corresponding Berry phase is non-zero 
and zero, respectively (cf., e.g., the two references ini&). 
This observation, too, indicates that the relationship be- 
tween the Berry phase and transport properties might be 
rather intricate. 

Before concluding, a few more comments are in order. 

In reality, the weak (anti)localization corrections are 
accompanied by the Aronov-Altshuler ones, which stem 
from the interference between disorder and Coulomb in- 
teractions. Despite a similar temperature dependence 
(5aa&xx — — {e 2 /irh) ln(7/T)), the former, but not the 
latter, gets quenched by an applied magnetic field, thus 
allowing one to discriminate between the two. 

Also, considering the potential attainability of the very 
low-density regime, it would be of interest to extend the 
analysis of the localization effects to the case \i < 7. 
However, as we have already pointed out, in this regime 
a systematic expansion of the conductivity corrections 
would only be possible for an (unphysically) large number 
of valleys. 

Moreover, even if this number were indeed large, the 
calculation itself would still pose a significant challenge. 
The main technical difficulty here is due to the double- 
pole structure of the Green functions (4), which results 
in the emergence of additional gapless Cooperon and dif- 
fusion modes in the RR (AA) channel, alongside the con- 
ventional (RA) one. A similar situation has been encoun- 
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tered in the context of dirty d-wave superconductorsiLiSi 
where it has been argued that the quasiparticle conduc- 
tivity receives aditional logarithmic contributions from 
the processes involving RR (AA) Cooperonsi^. 

Furthermore, the authors of RefiiS found the overall 
logarithmic correction to the spin (thermal) conductivity 
of normal ci-wave quasiparticles to be negative, contrary 
to (in this case, ill-defined due to a non-conservation 
of charge) electrical one. In view of the formal differ- 
ences between the Dirac-like descriptions of planar d- 
wave superconductors and graphene, it is, at least, con- 
ceivable that the correction to the electrical conductivity 
of graphene could change sign in the low-density regime, 
too. 

Lastly, we comment on the possiblity (as well as the 
need) of including other types of randomness. While the 
vertex (2) is devised as a general model of short-range 
(screened) impurities generating a random scalar poten- 
tial, it obviously misses out on those types of disorder 
that can be best represented by either a random vector 
potential or a random mass of the Dirac fermions. 

The latter might be relevant in the situation where a 
spatially inhomogeneous charge or spin density wave or- 
dering sets in. In turn, the former can be utilized to 
describe extended lattice defects (dislocations, disclina- 



tions, and cracks) in terms of a random magnetic fliD& 
In Ref£, it was assumed that the corresponding random 
vector potential (rather then the flux) is (^-correlated in 
space (< A q A- q >= const). Although this model pos- 
sesses a conformal invariance and, therefore, can be an- 
alyzed in quite some detai l 1 ?! 2 ? , a more adequate to the 
task would obviously be the model of a long-range corre- 
lated random vector potential (< A q A- q >oc 1/q 2 )- The 
latter is known to show a rather different behavior even 
in the high-energy (ballistic) regime, as demonstrated in 
Ref^i. 

A further investigation into these and related issues 
will be presented elsewhere. 

In summary, we carried out a comprehensive analysis 
of the first-order quantum correction to the conductiv- 
ity of doped graphene. We identified the conditions un- 
der which the conductivity correction becomes positive, 
negative, or zero. The earlier analysis of Rcf. 8 has been 
critically assessed, and the open problems have been out- 
lined. 

This research was supported by NSF under Grant 
DMR-0349881. The author acknowledges valuable com- 
munications with E. Abrahams, A. Geim, V. Gusynin, 
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